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' Abstract 

D ■ 

' We study a rotating Bose-Einstein Condensate in a strongly anharmonic trap (flat trap with a 

finite radius) in the framework of 2D Gross-Pitaevskii theory. We write the coupling constant for the 
interactions between the gas atoms as l/e^ and we are interested in the limit e — > (TF limit) with 
the angular velocity O depending on e. We derive rigorously the leading asymptotics of the ground 
^ , state energy and the density profile when Q. tends to infinity as a power of 1/e. If f2(e) = f2o/e a 

• "hole" (i.e., a region where the density becomes exponentially small as 1/e — > oo) develops for flo 

' above a certain critical value. If ^(e) 3> 1/e the hole essentially exhausts the container and a "giant 

, vortex" develops with the density concentrated in a thin layer at the boundary. While we do not 

\C ' analyse the detailed vortex structure we prove that rotational symmetry is broken in the ground state 

' for const. I loge| < n(e) < const. /e. 

^ ■ MSG: 35Q55,47J30,76M23. PACS: 03.75.Hh, 47.32.-y, 47.37.+q 
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In recent years much effort, both experimental and theoretical, has been put into the study of vortices in 
rotating Bose-Einstein condensates, see, e.g., jCD| . the review [FS| and the monograph [X] where extensive 
hsts of references can be found. Most of the theoretical research is carried out in the framework of Gross- 
, - Pitaevskii theory, whose status as an approximation of the quantum mechanical many-body problem was 

rN I established in [LSY| for the non-rotating case and in [LS| for rotating systems. On the mathematical 

. physics side an important topic has been the vortex structure in the strong coupling (Thomas-Fermi, 

TF) regime in harmonic traps when the rotational velocity is scaled with the coupling in such a way 
that the number of vortices remains finite |AD( IIMI) IIM2| . General results on symmetry breaking for 
sufficiently large interactions or rotational velocities and in traps of arbitrary shape were proved in the 
papers [SH [SI] that are not limited to the TF regime. 

Recently, attention has focused on rapidly rotating condensates where the number of vortices is much 
larger than unity (see, e.g., pCHSC, SCEMCj for experimental results). Much of this research has been 
for harmonic traps, e.g., [Mi ABB, ABm\ IABN2|. IWBPH IWGBPj . where "rapid rotation" means a 
velocity close to the limiting velocity beyond which the centrifugal forces destabilize the condensate, but 
anharmonic traps (mostly quartic plus harmonic) have also been discussed |AAB[ [El IBPi [0 IFBl IFZi 
IKBI IKFl IKTU| . For harmonic traps the eigenstates of a noninteracting rotating gas fall into Landau 
levels and rapid rotation implies, also for an interacting gas, that essentially only the lowest Landau level 
(LLL) is occupied. Using this fact detailed informations about the lattice of vortices have been obtained 
in [ABD| lABNl) IABN2|, IWBPl] . Some results for harmonic traps going beyond the LLL approximation 
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are discussed in jWBPl|, IWGBP] . In an anharmonic trap a restriction to the LLL is not adequate. The 
reason is that the energy gap between Landau levels is proportional to the angular velocity while the 
centrifugal energy is proportional to the angular velocity squared. In an anharmonic trap the latter can be 
much larger than the former, while in a harmonic trap close to the limiting angular velocity the potential 
and centrifugal energies almost cancel each other and the LLL energy is the dominating contribution. 

In the present paper we study a rapidly rotating gas in a trap that is as far from being harmonic as 
possible: The gas is confined within a finite radius R and the trap is "flat", i.e., the confining potential 
is constant (zero) inside the trap. Formally this trapping potential can be regarded as a limit of a 
homogeneous potential V{r) ^ {r/RY with s oo. Such a limit naturally leads to Dirichlet conditions 
at the boundary, but it is mathematically somewhat simpler to consider the case of Neumann (or free) 
boundary conditions and this is what we shall do. In this way the interplay between rotational effects 
and the nonlinear interaction terms are brought out in a particularly clean way. Dirichlet boundary 
conditions lead, in fact, to exactly the same results in the TF limit as we shall also show. Generalizations 
to homogeneous potentials with s < oo are in principle straightforward but the case s = oo merits a special 
treatment because it brings out clearly the essential differences between harmonic and anharmonic traps 
and also because of some special features with respect to the breaking of rotational symmetry. This will 
be discussed in Section 2.1. 

Our main results concern the density profile and the ground state energy in the asymptotic limit when 
the coupling constant l/e"^ (see below) tends to infinity (TF limit) and the rotational velocity il{e) is at 
the same time scaled with e. (The TF limit of the 2D GP functional without rotation is discussed in 
|LSY| .) Our estimates are not sharp enough to rigorously uncover the vortex structure of the condensate, 
but the variational functions that we use and which give the correct energy to leading order in e provide 
important hints about this structure. In particular, the regimes ^{s) 1/e and f2(e) S> 1/e require 
different variational functions, the former with a lattice of vortices distributed over the trap and the 
latter with a "giant vortex" in the region where the density is exponentially small. 

When considering the TF limit there is an important difference between traps that confine the gas 
strictly to a bounded region and traps where the gas can spread out indefinitely. If one considers for 
instance a trap given by an homogeneous potential V{r) = r**, for some < s < c» and performs the TF 
limit in a naive way, the result is trivial, namely the minimizer goes to zero and the energy to infinity. In 
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order to obtain a non-trivial limit it is necessary to rescale all lengths by e 2+= . In the case of infinitely 
high walls considered here the characteristic length of the problem is fixed from the outset and therefore 
no rescaling is needed in the TF limit. Consequences of this difference for the question of symmetry 
breaking in the rotating case will be discussed in Section 2.1. 

Rapidly rotating condensates in a flat trap have been previously studied by Fischer and Baym in 
[FBj and the paper of these authors triggered, in fact, the present investigation. Our analysis underpins 
and extends their general picture by rigorous estimates. We do not, however, confirm that the transition 
to the "giant vortex" state takes place for Q{e) ^ l/(e^|loge|) as imphed by Eq. (20) in [FBj . Our 
conclusion is rather that such a state emerges asymptotically at all rotational velocities ^^(e) ^ 1/e. The 
reasons for this difference are discussed in Section 2.4. 

We now define the setting more precisely. The starting point is the 2D Gross-Pitaevskii (GP) energy 
functional 



where Bb. denotes a ball (disc) of radius R centered at the origin, L the third component of the angular 
momentum (i.e. L — —id /d-d in polar coordinates (r, f?)), ri(e) the angular velocity and e is a nonnegative, 
small parameter. 




(1.1) 




(1.2) 




(1.3) 



By Sobolev immersion H^{Br) is contained in L^{Br), so the functional is well defined on f 



and denote by 'i'f^ a corresponding minimizeiH. Indeed, for any fl{s) < oo, one can prove (see for 
instance |Slj ) that the functional is bounded from below and there exists at least one minimizer. 

From the physics point of view a minimizer of (jl.ip describes the macroscopic wave function (the 
wave function of the condensate) of a Bose-Einstein condensate in the rotating reference frame. The 2D 
description is a simplification that is justified either in the limit of thin ("disc shaped") 3D traps, or traps 
that are very elongated along the rotational axis ( "cigar shaped" traps) so that the 3D wave function is 
essentially constant along this axis. In both cases the coupling is proportional to Na/h where a is 
the scattering length of the two-body potential (for its definition see, e.g., the Appendix in [LYj ) , N the 
particle number and h the extension of the 3D trap along the rotational axisH 

In the non-rotating case, ^{e) — 0, the minimizer is actually unique, by the strict convexity of the 
functional, and it is given by the (normalized) constant function, 
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The ground state energy is then 



If the angular velocity is different from zero the minimizer may not be unique since a rotational symmetry 
breaking occurs for f2(e) above a certain threshold value as will be shown in Section 2.1. 
We point out that the dependence on the radius R of the trap can be scaled out: 

i?fP = £;fP(i) = i?2£;GP(i?) (1.4) 

so that, without loss of generality, we can choose R = 1 and denote the functional by 
The GP functional can be rewritten in the following form that we are going to use: 

W ^ I {I (V - ,A,) *f ^ e(£&lffi + ffi} (1.5) 

where is the vector potential associated with the rotation, i.e., 

Mr) EE X f, (1.6) 

with z the unit vector in the ^-direction. In p.Sp one can recognize an analogy with the Ginzburg- 
Landau (GL) functional (see, e.g., [BR| ) in the theory of superconductivity. The vector potential (II. 6p 
is in this context due to a uniform magnetic field, while the wave function of the condensate is the GL 
order parameter (density of Cooper pairs). Using the L^— normalization of the minimizer, the analogy 
can be made even closer, namely the minimization problem in (II. 3p is equivalent to the minimization of 
the functional 
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over L2_]^Qi-]-Qa,lized functions. At this point, however, an important difference becomes evident, namely 
the presence of the centrifugal energy (the second term in the expression above) , which in the GL context 
could be interpreted as an electric field. This contribution, usually not present in the GL functional, 
is proportional to the square of the angular velocity and we are going to see that, in the regimes we 
are considering, it is responsible for a rather different behavior of the minimizer compared to GL theory. 
Another important difference between the GP and the GL minimization problems is the L2_jiQj-jnalization 
condition, that prevents, for instance, the minimizer from being identically zero, as it can be in the GL 
case. It also gives rise to an additional term (chemical potential) in the variational equation associated 
to (foil . 

In the next Section 2 we introduce some notations and state the main results of this paper. We first 
discuss the problem of spontaneous symmetry breaking in the ground state, then we study the regimes 
fl{e) < 1/e (Section 2.2), n{e) 1/e (Section 2.3) and fl{e) > 1/e (Section 2.4). Section 3 is devoted 
to the proofs, while in Section 4 we comment on the results and perspectives. 

^Any result for 'Sf^^ stated in the following is meant to be true for any minimizer, if it is not unique. 
^In "thin" traps a different formula applies at extreme dilution where the coupling becomes independent of a and depends 
logarithmically on the average density (see |SY) ). 



2 Main Results 



2.1 Spontaneous Symmetry Breaking in the Ground State 

The GP functional for a rotating 2D condensate in a general trap has already been studied in SI . A very 
interesting phenomenon generated by the rotation is the spontaneous breaking of rotational symmetry in 
the ground state. If the trap potential is polynomially bounded at infinity one can prove (see Theorem 4 
in [HI]) that for any fixed angular velocity there exists eq such that, if e < eo, no ground state of the 
GP functional is an eigenfunction of the angular momentum. The rotational symmetry of the functional 
is then spontaneously broken at the level of the ground state. An important consequence is that the 
minimizer is no longer unique, since a rotation by an arbitrary angle gives rise to a state with the same 
energy. 

A crucial ingredient of the proof in \Sl\ is that, in a polynomially bounded potential trap, the density 
of the minimizer tends to zero as £ ^ 0. In fact. Theorem 4 in [ST] is not true in the case of a trap with 
infinitely high walls as we are considering and we actually expect the opposite behavior: If fl is kept 
fixed, then for e sufficiently small the ground state is an eigenfunction of the angular momentum, and 
after an appropriate choice of a constant phase factor, a unique, strictly positive radial function. This 
difference can be understood by noting that in a trap of radius R the kinetic energy of a vortex is of the 
order R~^\ \oge\ for small e. Thus, if R is fixed, an angular velocity of order | loge| is needed in order to 
create vortices. In a polynomially bounded trap, on the other hand, the effective radius of the condensate 
increases as £ ^ and the critical velocity for the creation of a vortex behaves as e'^/(*+2)| loge| for a 
trap potential ~ r*, cf. the remark at the end of Section 3 in [ST]. Any fixed fl thus exceeds the critical 
velocity as £ — s- if s < oo. 

Despite this difference, our proof of symmetry breaking is obtained partly by a modification of the 
arguments of Theorem 2 in [Sl] . The following Proposition 12.11 states that for angular velocities smaller 
than l/y/ne, symmetric vortices of degree higher than 1 are unstableQ. 

Proposition 2.1 (Instability for Higher Vorticity) 

Let ^'„(r), n > 2, be the unique minimizer of£'^'^['^'] on the subspace of functions with angular momentum 
n, I.e., on any ^{e) < 1/ ^/tte, is unstable, i.e., it is not a local minimizer 

Proof: From the variational equation satisfied by the radial part of = f„(?')e*"'', 

-Ae„ + ^ - f^(£)<« + ^ = l^n{e)^n 

where the chemical potential /!«(£) is fixed by the L^-normalization of it is not hard to prove by 
a rearrangement argument (see, e.g.. Lemma 1 in SI ) that ^„ is a positive non-decreasing function, 
^nir) = 0{r'^) as r — > and ^^(1) = 0, i.e., £,n satisfies Neumann boundary conditions. Moreover by a 
subharmonicity argumenlH we can also prove the bound: 

< y {finis) + nie)n - n^} (2.1) 

Indeed, suppose that n > 1 and the opposite is true, then setting 

enir)>e' {fin{e) + n{e)n-n^)/2] 

we can have two possibilities: Either B'' ~ 0, and then the result easily follows, or it is an open interval, 
= (R^ , 1), by monotonicity of and A^„|b> > 0. In this case, by integrating A^„ over B^ and 
using Neuman boundary conditions, one has 

/ Ae„rdr = -i?>C(i?>) > 
Jb> 

*In the opposite regime of weak coupling and fixed Q, vortices of degree 2 or higher may be energetically favorable in 
anharmonic traps ILul . 

^For a similar proof see, e.g., Lemma 2.1 in [L5Y] . 



6> = |r e (0,1) 



which is a contradiction because ^„ in non-decreasing. 

The rest of the proof coincides with the proof of Theorem 2 in [ST]. Using the estimate (|2.ip . we 
can extract, as in (2.33) in [ST], a sufficient condition on the chemical potential for instability of the 
corresponding vortex: The symmetric vortex of degre^l n > rf G N is unstable if 

-{d - If fie + {d? - l)^{e)n -{d- \fn^ < 

or, choosing d = 2, 

-^„(e) + 3ri(e)n-n^ < 0. 
From the definition of the chemical potential and Schwarz's inequality it also follows that 

JBi i ^ ^ ) "^^ 

Inserting this bound in the condition above, we have instability if 
Hence any vortex of order n > 2 is unstable, provided ri(e) < l/^/ne. 

□ 

From Prop. [211] it follows that in order to prove the symmetry breaking in the ground state for a given 
^(s) l£ l/y/TTS it is sufficient to show that a rotationally symmetric vortex of degree smaller or equal to 
1 cannot be a minimizer of the GP functional at this angular velocity. This in turn can be achieved by 
exploiting some energy estimates. Let us first define 

Enie)= min / dr [ \V^\' + ^ + ^\ 

and 

so for fl{e) > Slft(£) no symmetric vortex of degree n < n can be a global minimizer of the GP functional. 
Since finis) < {2n + l)0o(£) (see Lemma 3 in [Sl]) we can use an upper bound on Ei{e) — Eo{e) to prove 
the symmetry breaking. 

Proposition 12.11 applies only for il{e) < I/^/tte but symmetry breaking can, in fact, be proved for 
fl{e) < C/e with an arbitrary constant C by using Theorem 12.11 that is proved later in the paper. Hence 
a part of the proof of the next proposition will be postponed to the end of Section 2.3. 

Proposition 2.2 (Symmetry Breaking in the Ground State) 

For e sufficiently small, no minimizer of ['^] is an eigenfunction of the angular momentum, if 

ellogel + 3 < n(e) < — 
e 

for any constant C G MJ^ . 

Proof: Using the normalized trial function 

{- if < r < e 
1 otherwise 



®The variational parameter d £ N is involved in the definition of a suitable trial function used in Theorem 2 in | S1| . The 
requirement n > d is necessary, otherwise such a trial function does not belong to H^{Bi). 



we can prove the upper bound 

i?.(e)<^+log(^)+l. 
Since ri„(e) < {2n + l)ilo(e) (see Lemma 3 in [ST]) and Eo{e) ~ we get 









< 3 




+ 1 



Hence no symmetric vortex of degree < 1 can be a global minimizer of the GP functional if > 
6 1 loge| + 3. On the other hand vortices of degree higher or equal to 2 are excluded by Proposition 12. II 
provided that ^{e) < l/^/ire. The proof of symmetry breaking for general il(£) ^ 1/e will be given at 
the end of Section 2.3. 

□ 



2.2 Energy and density for VL{e) ^ 1/e 



If <C 1/e, the rotation has no leading order effect in the TF regime. More precisely the energy 
asymptotics is the same as for a non-rotating condensate, and the density profile, l^f^^p, converges to 
the normalized constant function, namely the minimizer of the GP functional without rotation: 

Proposition 2.3 (Energy and Density Asymptotics) 

For any f2(e) such that lim£_>o efi(e) = and for e sufficiently small 

e^E^^ = 0{e^n{e)^) (2.2) 

TT 

Proof: Since ||5'^^||i2(gj) — 1 and i? = 1, we have 



and ||*?^|li4(Bj) > l/"" by Schwarz's inequality. Hence ()1.5|) leads to the lower bound, 

^2^GP > 

" ~ n 4 

The upper bound is obtained by evaluating the functional on the trial function namely 

e^RGP < i. 

TT 

Moreover, since ||5'^^||x,2(ei) = 1, this estimate implies 

iiiTrPi? , ii2 e'^rt(e)'^ 

||l*e I ^'^\\l^{Bi) - 4 

The L^— convergence of the density profile now follows by Schwarz's inequality. 

□ 

We stress that the result above says nothing about the fine structure of the minimizer and also nothing 
about its uniqueness. As far as the density profile is concerned, the first critical velocity at which some 
new effect comes into play is ri(e) ~ 1/e as it will be discussed in the next subsection. On the other 
hand, the fine structure of '^f^ depends on the angular velocity, even if il{e) <^ 1/e. If n{e) is simply a 
constant and e is sufficiently small, it is not hard to see that the minimizer is unique. More precisely, it 



is a radial function (and hence an eigenfunction of the angular momentum) which can be chosen strictly 
positive. In this case the result in (|2.3p can be improved and the convergence can be extended to L°°{Bi). 

According to the discussion in [AD] and the rigorous analysis in [IM1[ IIM2j of rotating Bose-Einstein 
condensates in harmonic traps, the first critical velocity for the occurrence of vortices, i.e. isolated zeros 
of the minimizer, is in that cas^ of the order ^.{e) ^ e| \oge\. More precisely, if Vtd{e) < il{e) < ^l,d+ii£), 
where 

n^ie) = ce[\ logel + {d-l) log | log ej] , (2.4) 

the minimizer has exactly d vortices of degree 1. A similar behavior was shown in [Se' for a slightly 
different model of superfluids. 

Such results together with the considerations in Section 2.1 suggest that in a flat trap vortices start to 
occur if r2(e) ~ I loge| and the rotational symmetry can be broken. The spontaneous symmetry breaking 
cannot be seen at the level of the density profile IvP^^P however, because the average size of each vortex is 
very small (area of the core of order e) in the TF limit. The total vorticity of the minimizer is proportional 
to the angular velocity, provided that ri(e) ^ |loge|, and therefore, as long as r2(e) <C 1/e, the region 
covered by the vortex cores has Lebesgue measure zero in the limit e — > 0, in accord with (|2.3p . 



2.3 The Regime Q{6) ~ 1/e 

In the regime il{e) ^ 1/e the rotation is so fast that it modifies the density profile itself: Since the 
centrifugal energy in (II. 5|) is of the same order of the non-linear term, it is no longer convenient for the 
condensate to be uniformly distributed over the trap, like in the non-rotating case. Such an effect can 
be seen at a macroscopic level, namely the density profile converges to a non-constant function, which 
minimizes a TF-like functional. 

Before stating the main results we first need some new notations. Without loss of generality we can 
assume that f^o = eil{e) is a constant independent of e. Moreover, for any fio > 0, we introduce the TF 
functional, 



df 



defined on the domain 



V^^ ^{peL\B,)\p>Q]. 
The functional above has a unique minimizer, p"^^ , and we denote 

pTF _ 

pev 



mm 



(2.5) 
(2.6) 

(2.7) 



The minimizer p can be explicitly calculated: 



1 

TT 
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^(l-2r2) 



a 



no 



if no < 



if no > 



(2.8) 



where [ • ]+ stands for the positive part, and the ground state energy is 
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E 



iTF 



1 

TT 

no 

4 



"0 



nnf) 

768 

-r!o 



if no < 



if no > 



(2.9) 



^The overall factor e in the critical velocities 112.41 1 is due to the scaling mentioned in Section 1 and Section 2.1. 



We point out that, if Qq > p^^ has a "hole", i.e., a macroscopic region where it is identically 
zero, centered at the origin: With 

we have p^^{r) = 0, for any r < Rq. We also define 

Po = supp(pTF) = {re6i|r>i?o}. (2.11) 
The first result concerns the energy asymptotics: 

Theorem 2.1 (Energy Asymptotics) 

For any fig > and for e sufficiently small 

^2 ^GP ^ ^TF ^ jQg^i)^ (2.12) 

The leading order term, proportional to 1/e^, in the asymptotic expansion of E^^ is due to the 
centrifugal bending of the profile while the remainder (of the order | loge|/e) is the contribution coming 
from the fine structure of the minimizer. Indeed, is expected to carry a very large number (of the 
same order as ^{s)) of vortices of degree 1. As suggested by the trial function (|3.2p used in the proof of 
Theorem l2.1l fsee also |FB| ). such vortices should be distributed over a lattice with a spacing of order ^/e, 
so that the average vortex core covers an area proportional to e. A simple argument (see, for instance, 
[BBH2j ) shows that the kinetic energy of each vortex is of the order | loge|. This explains why the total 
energy contribution of vortices produces a remainder of the order | loge|/e. 

As stated in the Introduction, the results proved in Theorem 12.11 and in the rest of this Section also 
hold in the case of Dirichlet boundary conditions (see the Remark |3. II in Section 3.1). The crucial point 
is that the limiting functional (|2.5p contains no kinetic energy and hence boundary conditions become 
irrelevant in the TF limit, at least to the leading order. 

The convergence of the profile j^f^^p to p^^ is a straightforward consequence of Theorem 12. II 

Corollary 2.1 (Density Asymptotics) 

For any flo > and for e sufficiently small, 



- P^niLi(B.) = OiV^lT^l). (2.13) 

If JIq > the estimate above can be improved and we can prove that the profile j^f^^p is exponen- 
tially small in e inside the "hole" , i.e. where p^^ is zero: 

Proposition 2.4 (Exponential Smallness of the Density in the "Hole") 

Denote 

r, EE jj^e 6i I r < i?o -e^} (2.14) 

where Ro is defined in (|2.10p . For any flo > and e sufficiently small, there exist two constants 
and C'q^ such that, for r ^T^, 



< CoosVllogel exp 



C'^Aisi{r,d%f 



£3 

The results stated above allow us to complete the proof of Proposition 



(2.15) 



Proof of Proposition 12.21 

It remains to prove the statement for any 51 (e) = ^o/s- Suppose that the opposite statement is true, 
namely the ground state energy is reached on a symmetric vortex of the form ^^6™''- Then there 
must be some £ N such that 

E,Ae) - = + (2.16) 



where we have used the upper bound for E^^ proved in Theorem 12.11 Using the rough lower bound 
Enis) > n^, we immediately get the upper bound < C/e for some constant C independent of s. 
The right hand side of (|2.16p can be bounded below by 



/ df 


n fi(e)r 




r 2 



<7TFrc2i 



[ df 


n 




2 



n 2 



E 



iTF 



Therefore one has the estimate 



df 



C\\oge\ 



Since ^n^e™^'' is a ground state, it must satisfy the estimate (|2.13p . and then 



C|l0g£| 



> / df 

Bi 



n 2 



etM> / df 

'Bi 



Cv^logi 



£2 



or 



df 



where we have used the bound 



df 



n{ey 



n 2 



n{ey 



Coo\/|log£| 

3 
£2 



(2.17) 



TFI 



> 



£2 



The left hand side of (|2.17|1 can be explicitly calculated: If flo < 4/^/7r, one has 



df 

while, if r^o > 
df 



12 



n{ey 



By minimizing over rig, i.e., taking 



in the first case and 



in the second, we get 



_ m 
4 



1 + 



df 



Bi 



n{ey 



p'Hr)> 



72 72£2 



if ^^0 < 4/^71^, and 



df 



Bi 



n{ey 



p''{r)> 



1 32 



3 971^1 



£ 



2 ■ 



if Hq > ^l\fTi- Therefore in both cases ()2.17p implies that 



0<C^,„ <Co„v/^logi 

for some strictly positive constant G^^^. For £ sufficiently small this is a contradiction and then no 
symmetric vortex can be a ground state of the GP functional. 



□ 



2.4 The Regime (^(e) > 1/e 

In order to present the results in a transparent way we assume that the angular velocity is a power of 
namely 

for some a > cH. In this case the limiting functional is analogous to the one introduced in Section 2.3, 
provided fio is replaced by Oi/e" and the energy scale by e^", i.e. 



f-'w-e-/ dnp'~%^\- (2.18) 



The ground state energy of this functional, i.e. 



min 5r[p]=£r[pr] (2.19) 



is given by 

and the corresponding minimizer is 
where 



02 / Sc-CK 

rpp ^ "1/1 



pr = S^[^^-i??], (2.21) 



8e' 



Hence the function p^^ is supported in a very thin layer near the boundary and, as e 0, it converges 
as a distribution to a radial delta function supported at r = 1: If F(r) is a continuous function, one has 



dTpY{r)F{r) = ^ C ""'^ dzzF(^7+m) = f\z z F { J + rA ^ F{1) 



8£2"7o ^ Jo V V 7'^" 

We can now state the main results for this regime, starting with the energy asymptotics for e ^ 0: 

Theorem 2.2 (Energy Asymptotics) 

For any fii > 0, a > 0, and for e sufficiently small 

^2+2a ^GP ^ ^TF ^ Q(^^2a^ ^ ©(e^j logej). (2.23) 

The two remainders in the asymptotic estimation of the GP energy have different sources: The 
second one, of the order | logej/e^", is due to the convergence of the density to a delta function and the 
radial kinetic energy that is ignored in the TF functional. The other term, of order is due to the 

approximation of the vortex structure in the region where the density is exponentially small by a trial 
function with a single vortex located at the origin. It is clear that the estimate ()2.12|) is not the a ^ 
limit of (|2.23p . We also note that for a > 2 the last error term in (I2.23P is larger than the second term 
in (H^. 

Since the function pj^ does not converge in any L^'-space, a result analogous to CoroUarv 12.11 does 
not hold. We are able to show, however, that the L^-norm of '^f^ converges to zero almost everywhere, 
except for a thin region (with a size of order of a suitable power of e) near the boundary. 



'Our analysis applies, in fact, to arbitrary angular velocites 0{e) 3> 1/e, one just has to replace Qi/e°' by ef2(e). 



Corollary 2.2 (Density Asymptotics) If > and e is sufficiently small, 



= 0(£") + 0(£2-"|log£|) 



for < a < 2, while for a >2, 



(2.24) 
(2.25) 



with R,^i3 = (1 - e'^y^'^, for any 1 < (3 < 2. 



The above estimate is strengthened in the following proposition. The reason why we state Corollary 
^separately is the analogy with the previous Corollarv l2.1l It is also used in the proof of the following. 



(2.26) 
(2.27) 



Proposition 2.5 (Exponential Smallness of tlie Density) 

Denote 

r; = |r e Bi\r< 1 -e"'/^} 



where a' = min [a, 2 — a] and (3 is any number such that 1 < (3 < 2. 

For any fii > there exist constants Cui and Cq_^ , such that for e sufficiently small, 



|«'fP(f)|' < Coie"'/'|log£| exp 
if Q < a < 2 and f ^ , and 



C^^dist(f,9T;)2 



exp 



c^^dist(f,ar;')^ 



(2.28) 



(2.29) 



if a>2 and re T^' . 



A straightforward consequence of the above estimates together with the normalization of ^'^^ is that 
the density of any minimizer of the GP functional converges to (5(1 — r) in a distributional sense, in accord 
with the discussion in [FB]. 

Another important difference compared to the regime f2(e) ^ 1/e is the form of the trial function 
p.36p used in the proof of Theorem 12.21 This function is an cigcnfunction of the angular momentum, 
i.e., the whole vorticity is concentrated at the origin. On the other hand, Proposition 12.21 implies that 
the true minimizer cannot be an eigenfunction of the angular momentum, at least as long as fi(e) ^ 1/e. 
Nevertheless we expect that the number of vortices contained in the region where ^'^^ is not exponentially 
small is negligible compared to the total vorticity of the function (see, e.g., the numerical simulations 
contained in ^KTUj)- A wave function of this kind is often referred to in the physics literature as a "giant 
vortex" . Since the vortex contribution to the energy depends essentially only on the winding number at 
the boundary of the thin region where the wave function of the condensate is not exponentially small, 
a trial function with the vorticity concentrated at the origin can lead to a good approximation to the 
energy. 

This behavior is also suggested by the fact that the minimization of a modified GP functional over the 
subspace of functions with fixed angular momentum with a subsequent minimization over the value of 
the angular momentum gives a ground state energy with the same leading order asymptotics as Ej^ and 
Ef^. Indeed, if we define (c.f. [FB] ) 



df 



and 



E. 



.tf' 



Bi 



= mm 



min 

P e v^^' 
Jp=i 



cTF' 



[P] 



(2.30) 



(2.31) 



where T>'^^' is the natural domain for the functional (|2.30p . then it is not hard to see that 



hm E^' = hm E^^ = lim e^+^^ijGP ^ (2.32) 

even though Ej^' > for any e > 0. The functional (|2.30p is obtained from p.Sp by neglecting 

the radial part of the kinetic energy and restricting it to eigenfunctions with fixed angular momenturqfl 
^l^i+a ^ In other words the TF functional (|2.30p describes the asymptotic behavior of the GP functional 
almost as well as (|2.18p . 

However, while the heuristic discussion in |FB] suggests that a giant vortex occurs only for angular 
velocities larger than l/(e^| logej), we found no evidence in our rigorous analysis that a change in the 
minimizer occurs above that threshold. In fact such a critical angular velocity for the transition to the 
giant vortex is estimated in jFBj by imposing the condition that the ground state energy Ej^ equals, 
to leading order, a certain upper bound for This upper bound, however, is calculated in [FB| 

with a trial function of the form (|3.2p and is not optimal. In fact, the same comparison with a better 
upper bouncf^ for e^+^"i?^^ gives the correct answer, namely that £^+^"i?^^ is close to Ej^' for any 
a > (see (|2.32p ). Hence the transition to the giant vortex should occur for any angular velocity 51 (e) 
of order higher than 1/e. 



3 Proofs 



3.1 The Regime n{e) ~ 1/e 

The main result concerning the regime ~ 1/e is Theorem 12.11 and we start by proving it. Some 
technical but crucial details of the proof are contained in Subsection l3.1.1l f Proposition 13 . 1 1 and Theorem 
13. ip , where we present some estimates for the kinetic energy of the trial function. 

Proof of Theorem 12.11 

We are going to prove the result by comparing an upper bound for the ground state energy with a suitable 
lower bound. 

Lower Bound: The lower bound for E^^ is actually trivial. By simply neglecting the positive contri- 
bution of the magnetic kinetic energy in (jl.5p we immediately get 



TF 



(3.1) 



Upper Bound: We prove the upper bound by testing the functional on a trial function of the following 
form 

^[r)^c,U{r)xe{7^9e{r)- (3.2) 

The radial part is given by 



„TF 



if VIq < 



fe{r) 



(3.3) 



3e^/p^ if flo>^ 



where je is a suitable cut-off function to regularize \Jp^ at the boundary of the hole. Our choice is 

r if r < i?o 



1 



if i?o < < -Ro + £ 



otherwise. 



(3.4) 



®More precisely, the correct value of the angular momentum should be the integer part of i^/e^+" but the difference 
between these two quantities produces a correction of smaller order in 112.311 1. 
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For instance the one calculated using the trial function H3.36| l in the proof of Theorem 12.21 



The function is a phase factor that can be expressed in complex coordinates z = x -\- iy as 



=w = n 



\Z - Zi\ 



(3.5) 



where £ is a square lattice of spacing defined in the following way: 



< 1 - 2V2ie 



(3.6) 



We assume that the spacing is of order i.e. £^ — for some 6 > independent of e, so that the 
number of lattice points, denoted by A^g, is proportional to 1/e. 

We note that the phase carries vortices of degree 1 centered at the lattice points. Moreover, each 
vortex core is contained inside the fundamental cell and its radius is of order y/e. This choice is suggested 
by previous works (see e.g. |BBH2| liMl) IIM2] ) on rotating condensates, where it is shown that vortices 
of higher degree than 1 are energetically unfavorable. 

Since g^ is not differentiable at the points of the lattice we need to multiply it by a function, Xe, that 
vanishes at these points and we take 



Xe{r) 



if |f-fJ>e'' 



if Ir-fJ <£'' 



(3.7) 



for some rj > 1/2. 

Finally the constant is fixed by the normalization condition and it can be easily checked that, for e 
sufficiently small, 

1 < c? < 1 + Ce. 



Setting 



A = Si\ U 



where i3* is a ball of radius centered at f^, the functional evaluated on the trial function 

by 



(3.8) 

(3.9) 
is given 



J A J A 



e2 



dr 



2 



ci / iV/.r + c^ / drf^ 

J A J A 



+ ct I dr\W{Xefe)\ + 



+Ct I dfxlfe 



<ci I dr |V/, 
/8i 



(V - iX) 9e 



2 gTFu^ 



< 



drxlfe 



<ci / df\Vf,\' + Ci / df 

JBi J A 

+ C3 



(v - iA,) 







[ dr 


Ae 







\/~iA,] g, 

2 £TF[|^|2] 



C2 / drxl iVffe 



— < 



Cj C 

^ £2 + ^3-4„ + £ 

where we have used the uniform boundedness of fe, the estimate p.Sp . and the fact that the number of 
lattice points iVg is bounded by C/e. 

The gradient of the phase can be bounded from above inside any baU Bl: 



< 



1 



|Vg I < V 1 < 1 ^ 

~ \r-rl\ Mj^i\r-fj\ " \f-fl\ ' £^ 



for any rG B^, so that 

/ drxl iVffel' < 



(3.10) 



and hence 



< / df|V/e| +Ci / dr 



[ dr 




-) 5e 


Ik 







2 £:TFr|^ 



for any 77 > 3/2. 

Moreover, the radial part of the kinetic energy can be bounded by a constant if VIq < and by 



d^|V/,r<Ci/ dr\\/j,\' + C2l drr4^<C3 + a 



Bl 



Ro 



dr 



<C|l0g£| 



if ilo > Then we get 



<Ci I df 



[ df 


{y - iA, 


-) 9e 


Ik 







£ 



for a possibly different constant C2. 

The upper bouncf^ now follows using Proposition 13.11 and Theorem 13.11 in the next section, choosing 



^ ^ \ W ^'^^ < 77 < 00: 



(3.11) 

□ 

Remark 3.1 (Dirichlet Problem) 

The proof in the case of Dirichlet boundary conditions, i.e., if (|1.2p is replaced with iJp (Si), looks exactly 
the same: The trial function has simply to be multiplied by a cut-off function, which is 1 everywhere 
except for a very thin region in the neighborhood of the boundary where it goes to 0, in order to satisfy 
the required boundary conditions. The error coming from such a cut-off function can then be included 
in the remainder in (|3.1ip . 



Proof of Corollary [211 

Let us first consider the case ilo < Using the explicit form of the TF minimizer p^^ (see ()2.8p ) and 
the estimate ()3.1ip . one can calculate 




Note that the constant actually depends linearly on r) (see the proof of Lcmma l3.1l in particular Eq. H3.32|l '). 



< 



and then the L^— bound foUows from Schwarz's inequahty. 
On the other hand, if iln > -?=, one has 



•Do 



8i\X)o 



+ f / drir' - Rl)\^T? < f - < Cn^M \oge\ 

where we have used again (12. 8p and p. lip . From the same inequahty one also has 

df(r2 -i?2)|^GP|2 < 



iTF I ' 



/ df (|M/,C=P|2_pTF)2^ /■ d^-'|,I,GP|4^^TF[|^GP|2]_^' 
Jvo JBi\Vo 

<£'^^[\^rf]-E^^<Cnoe\\oge 
f dr\^f^\^<Caoe\\oge\. 

JBi\Vo 



Bi\Vo 



SO that 



(3.12) 

□ 



Proof of Proposition [2741 

The proof is similar to the one of Proposition 2.5 in [AAB] . 
The variational equation satisfied by '^f^ is 



GP 



e £^ 

where the chemical potential /i^ is fixed by the L2_jiQi-]-naii2ation of '^f^: 



GP 



Setting Us = j^'j' | and using the simple estimate 



no 

e 



one can easily check that 



XI,GP*^^GP 



^AC/e < 



I^^GP|2 ^^§r-2|*GP|2 



t2 ■ 



Moreover, thanks to Theorem 12. II and Corollary 12.11 



„TF||2 



+ Cn„ Ve|loge| 



so that 



If we define 



2 



-2C/e + CooVe|loge| 



£3 



(3.13) 



(3.14) 



then, for e sufficiently small, the function is subharmonic in 7^ and therefore, for any point f G 71 
with dist(r, 97^) > g, 

C\\Us\\L^Bl\Vo) 



< 



Hence, using the estimate p.l2p and choosing, for instance, g = /2, we can conclude that 

Ueir) < CnoeVllogel 



for any r G 7^ . 
Let us now define 



For any f e 7^ , 



and 



i.e., U' is a subsolution in 7^ of 



CooE^vllogel 



.Ac/; + ^<o 

£3 



o<ui<i 



-Am + — = 

£3 



u{d%) = 1. 

On the other hand it is not so hard to verify (see e.g. Lemma 2 in [BBHlj ) that the function 



exp 



is a supersohition for the same problem if 



4£3 (i?o - £3) 



and hence for any e sufficiently small. The result now follows from the comparison principle. 



□ 



3.1.1 Technical Estimates 

In this section we want to present some estimates involving the function (|3.2p . We start by stating a 
simple but important result: 

Proposition 3.1 (Upper Bound on the TF Energy) 

Let 'J be the function defined in p.2p . with rj > 1 and JIq > 0, then, for e sufficiently small, 



£^^m']<E^^ + Cnoe. (3.15) 
Proof: A simple estimate using p.8|) shows that 



^TF[|^|2]„^TF[^TF]< J^dr \ [ct f^ - (p^^ 



4 



+ C£2'^-l<C'||c?/|-pT^||^,(^^^+C£2''-l 

where A is the region defined in (|3.9p and we have used the fact that p^^ and are both uniformly 

bounded and the area |Uig£K^| is bounded by Ce^^~^. 

Now 

ll-e/e - P^1,HB., < J/r \p^^ - Clfl\ + C£^-l 



and, setting 



1 < Ce, the first term on the right-hand side is bounded by 



if fin < -4=. On the other hand, if fin > -4=, 

/ dr |pTF _ < 4 /■ d^pTF ^ ^^^^ 

where the last term is due to the cut-ofF function and Vq is defined in (j2.1ip . 
The main result contained in this Section is the following 



□ 



Theorem 3.1 (Upper Bound on the Vortex Contribution) 

Let be the function defined in (j3.5p . — 6\/e and flo > 0- There exists a constant Cno,s independent 
of e such that for e sufficiently small and 5/2 < < oo 



< 



2e2 



(52 



Cao,<5|log£| 



(3.16) 



where A ( depending on ) is defined in 



Remark 3.2 (Vortex Lattice) 

As far as the leading order of the GP energy is concerned, the vortex structure of the minimizer 4"^^ is 
not so important: The choice of a regular square lattice in (|3.6p is just the simplest for computational 
purposes but the result in Theorem l3.1l is expected to hold for any trial function with vortices on a regular 
lattice, provided that 6^ in p.l6p is replaced with the volume of the rescaled fundamental cell, which is 
the relevant parameter in the estimate. 



Proof: Expanding the expression in (|3.16p . we get 

2 



dr 



\/-iA,] g, 



dr \Vg,\^ + — I drg* {fx Vg,) 



(3.17) 



+ 



4e2 



dr^r^\g,\' 



The last term can be easily bounded from above by 



02 

4e2 



8i 



2 ^^0 



Using the fact that g^ = e*"^, where 



f^) = arctan 



y-Vi 



(3.18) 



the second term can be explicitly calculated: By applying Stokes's theorem, 

J drg* (r X Vge) = — — J (irrxV(?i = — ^ J dfV x (^r^Vci 



no 

' 2e 



/A 
„2\ 



ds ■ r^'^cj) = 



(3.19) 



Since for any r £ dBl, 



12 21 

\r - rA 



the last term in the expression above can easily be bounded by 



where we have used the estimate p.lOp 



inf, 



< 



ec \r - n 



(3.20) 



Since the lattice spacing is chosen to be equal to (Jy^, the number of lattice points satisfies the bound 



< 



and then 



2£ J dBl 



< 



p3— Jj 



(3.21) 



(3.22) 



Moreover the sum appearing in (|3.19p can be replaced by the integral over Bi: let and be the 
fundamental cell centered at the origin and at respectively, 



1 



1 



n — / drr — -TT drr = 



6 



so that, setting Ae = Bi \ (U^e^Q*), 



^1 JUiecQl 



N^e 1 

r- H — < — 

6 -ijj^ 



dfr^ 



1 



< 



TT (l-C"4)2(^-iV,£2) 



2e 



df + C < 



C 



(3.23) 



because the lattice is chosen in such a way that, for any i e £, < 1 - 2^/2^e• 
From inequalities p.22p and (|3.23p we then get (for any 77 > 5/2) 



df (r X V.9e) < 



2ef2 

but the number of points in the lattice can be estimated below as 



7r(l - 2V24)^ 



< 



C 



(see for instance Theorem 7.7.16 in [H]) so that 



£ 



dr (r x Vg^) < — —:r + ^ < 



2eR 



2(52e2 



(3.24) 



(3.25) 



The first term in (I3.17P is the most difhcult to estimate and we deal with it in the following Lemma [5TT1 
Altogether the three upper bounds then give the result for a possibly different constant Coq.s- 



□ 



Lemma 3.1 (Kinetic Energy of Vortices) 

Let Qe be the function defined in p.Sp . £e — S^/e and 77 > |. There exists a constant Cs independent of 
e such that for e sufficiently small 

f IY7 |2 / , Cs\loge\ , . 

Proof: Wc first notice the useful fact tliat 

df \Vg,f ^ [ dr ^ [ df V4>^ 

Ja Ja 

where </> is defined in (|3.18p and cj) is the function 

4>{r^^J2^n\f-n\. (3.27) 

iec 

Indeed, 4> and are conjugate harmonic functions (the real and imaginary parts of In Yiii^ ~ ^i))) that 
doccj) = -dy(j), dy4> = dx(t>- 

Since 4> is harmonic, the last integral can be explicitly evaluated by means of partial integration: 



dr 

A 



2 



dBi on f^JdBi on 



where n stands for the outer normal to integration path. 
We are going to consider the two terms separately. 

Outer boundary: The contribution at the outer boundary is given by 



where we have used the complex coordinate notation, z — x + iy. 
The first step in the proof is the replacement of the sum over i with an integral over 



v/2 

1? 



^ln|e*''-z,f - / dz \n\e'^ - ^ f dz In ^' 



Thanks to the choice of the lattice p.6p , 

|e*''-z, -z| 1 



|2 



> 



le'^'-z.l 2 

because, for any •& € [0, 27r], \e^^ — Zi\ > 2\/2i^ and \z\ < £^/^/2. Using therefore the bound 

ln(l + i)>t-t^ 

which holds true for any t > —1/2, we get 

^ittc^Q^ ' |e^^-z,|' " ^^I^^S. I |(cosz9,sin,9)-r/ 

[-2r • ((cos I?, sin t^) - Fi) + " ^ 



|(cos?9,sini9) — r^l'* 



< 



< ^ V / ^^4[f • ((cosi?,sini9) - f;)]Vr4 ^ 
" ^ 1^ is. |(cos^?,sin^9)-f;|^ - 



iec 



4^2 ^4 

I (cos??, sin??) — fi\^ I (cos??, sin??) — fi\^ 



< 



I (cos??, sin??) — ■fi'f' |(cos??, sin??) — ri\'^ 



Since the functions 1 /r^ and are positive and subharmonic we can easily bound the expression above 

by 



I (cos ??, sin — Tif I (cos ??, sin ??) — ri j** 



< 



< 



dr 



+ 



(cos??, sint?) — |(cos??, sin??) — r|'' 



< 



< 



dr 



✓2 



(cos??, sin??) — |(cos??, sin??) — r|'' 



< C 



so that 

On the other hand 
1 



2:7; < 



dz In 



C. 



dz In 



4' 



dz In 



e — z — 



dz In 



V2 



dz In e**' - z < 



C|ln4| Si_M.\U,6£Q^ 



V2 



C|ln4 



< 



'(i-tr 



< 



qin4 



< 



where we have used the estimate p.24p for the number of points and the fact that 



dz ln|e^''-z| =0 



Br 



for any < i? < 1. 
Since the function 



a(z) 



2 - ze 



is positive for any ?? G [0, 27r] and z G '^i-se/v^i ^^'^ initial expression in (j3.28p is bounded from above by 



dBi 



dg.^^<iy d?? 



7^e(^?) 



where 



dz In 



Ie'^-z| 



V2 



is easily proved to satisfy the upper bound 

C|ln4| 6i_M.\U,e£Q^ 



V2 



< 



qin4 



We need now to replace the sum over j with an integral over B-i_2^^^ ■ Since the function a{z) is harmonic, 
we can apply the mean value theorem to get 

2 

For any j e £, the right hand side can be easily estimated using Harnack's inequality: 



Qe\Be 



1 



c 



< 



< C£ea{zj). 



In the same way it is possible to show that for s sufficiently small, there exists a possibly different constant 
C such that 

be{Zj) > -C£ea{Zj) 

so that 



and then 



Since 



< CL 



/ dz a{zj + z) = dz a 



and 



we conclude that 



< / rf^ a(2) < 



27r 



E 



2 - zje-^^ - z*e''^ 1 

« _ ^.12 72 



e'" — z 



^2; 



2 - ze"*'' - z*e*'' 



l-2x/2<e 



< C4A^e. (3.29) 



On the other hand, using again the harmonicity of a, one has 



dz 



2 - ze-'^ - z*e'^ _ 27r (l - 2\/24)^ 



Since for the fundamental solution of the Laplace equation 



1 

2^ 



d-& ln\Re''^-x 



InR if |f I < i? 
In Ixl if Ixl > R 



we then get altogether 



ds- ^ < - 
on 



27r(l-2\/24)' 



Jb^ 3i^\UjQi Jo 
V5 



(3.30) 



C 



dz I dd In \e''* -z\ + I d^ |7^e(^?)| + < 



<2. .up |«U<,)||«.«,)| + fZ<3|M,^i<Cd|ld 

i?e[0,27r] ^£ ^£ £ 

Inner boundary: A straightforward calculation gives 



(3.31) 



-V / ds ■ ^4> = 2TTr]NJlne\- y f dd Inle'^e'^ + Zi - zA + 



< 2T:r]Ns\\ns\ - 27r ^ In | 



< 



< 27r77iVe|lne| - 27r ^ In | 



, , Cem^\loge\ ^ Cs\\oge\ 



(3.32) 



where we have used p.20p . 

In order to get the desired estimate we need now to replace the sum over one of the two indices in the 
expression above with the integration on a suitable domain. Therefore the quantity which has to be 
estimated is the difference 



(3.33) 



Using the estimate ln(i) < ^(t^ — 1), which holds for any t > 0, we can bound the expression under the 
sum in the following way 



r,- — + r — r ' 



- 1 



E 



2it\n~rj\ JQ, 



df / dr' 



\f~ r'r + 2(r-r') • (ri - r,) 



y i J / dr f d/ (r'+r'^) < V 



where we have used the central symmetry of the fundamental cell Qe and the lattice C. 

On the other hand, since the function l/r^ is subharmonic and positive, one can easily prove that 



< 



c 



1 



r," — r, + r — r' 



< 



dr / dr' 



Si "'S'e |r-r' 



|2 - 



C r C^lloge 



< 



so that the difference in (|3.33p is bounded by Cs \ loge\/e. 

In order to extend the integration to the whole disc Bi , we observe that 



-y 



Qi JQl 



dr i df' In Ir — r^l < 



dr 



df' In Ir — 7^ 



[ df [ df \n\f-f\ + [ df j df\n\f-f\-^j df [ d7^1n|f-?='| (3.34) 



where stands for the domain Bi \ ^iecQ\- 

The last term in the expression above is bounded by 



1 



e J A. 



— I df I df \D.\f-f\ < I df I df* Inlr-r'l < - 



Ae 



e J A 



Cln4 ^ Cillog^l 



where Ae = Bi\ fi_2V 24- 

For the second term in (|3.34p . we can use p.30p to get 



f dr f df' hi|r — r*! = ^rr / df \ i dr'r' Inr 



Therefore one has from (|3.34p 



dr'r' Inr' )■ < 



^E/ drj dr^\n\r~r^\<^^J dr ( 



rfr m I r — r I H < 



< 



df [1- r'^\ + ^^I^Qg^l < ^ + ^^I^Qg^l 



so that 



and finally 



o , I I ^ , Cs\\oge 
-2tt 2^ In |z, - l < ^ + 



El 7 1^ 



, Cs\loge\ 



dBl 



(3.35) 



Combining this result with the estimate for the contribution at the outer boundary, we complete the 
proof. 



□ 



3.2 The Regime Vt{e) > l/e 
Proof of Theorem [2721 

The lower bound can be proved in the same way as in the proof of Theorem 12.11 so that one easily gets 



TF 



*(^) = Ceje{r) \ pj^ir) exp < i 



For the upper bound we evaluate the functional on the following trial function 

where we used polar coordinates, r = (r, t?), [ • ] stands for the integer part, is the cut-off function 

r if r < R, 

r^-Rl 



(3.36) 



ie{r) = < 



if i?? < r2 < i?2 + £/3 



(3.37) 



1 otherwise 

with some (3 > a, and is a normalization constant satisfying the following bounds 

1< 2^ < l + C£2^-2". (3.38) 

A simple calculation shows that 



+ 



(3.39) 



The first term in (|3.39p is bounded by (using (|3.38p ) 









[ df 






hi 




Jb 



dr 



< 



t2P 



Re 



2p 



,TF 



,TF 



dr 



< 



+C3e 



2a 



dr— 
~ r-^ 



drr{r'^ - Rf) 

2 



dp 



,TF 



dr 



dp, 



,TF 



m \ dr 



< 



-2q ^-2(1 



(The constants depend on the choice of (3 and C3 — > cx) if ^ 00.) Moreover, using p.38p and the fact 
that 



for some < < 1, we can estimate the second term in p.39p as follows 
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In a similar way one can prove that 
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drr3(l-,;2)pJF(r)< 
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and then the result follows if we choose a finite /3 > 2a. 



□ 



Proof of Corollary [272] 

Let us start by considering the case < a < 2. Defining T)^ — Bi\ Bu^ we first notice that for any non 
negative function p e L'^{'D^), normalized to 1 in L^{'D^), 

er [p,T^e]>Er 

where Sj^ [p, V^] denotes the functional 



£r[p,'De]^e^''J^ dr[p'-^ 



Hence, setting p^ = l^'^^P and 

~ _ Pe 

Pe = n n 

WPellL^V^) 



we get 



>Er\\Pe\\mv.)-'--e 
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1 - 
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{\\P^\\l^{V, 



- 1 



1 

ll^'ellLi(X),) 
P£llLi(Pe) 



> 



-Wp^l^v,)^ ^ ^ 1^ 

where in the last step we have used Schwarz's inequality and the fact that ||pe||2,i(p ) 
On the other hand from the upper bound in the proof of Theoreni l2.2l one has 

[Pe] < °P < Er + Cie'" + €'26' \ loge\ 

and then (omitting for simplicity the subscript L^(I?e)), 

Er \\Pe\\ + (1 - llPell) < Er + Ci.^" + C^S^I log 

and therefore 



< 1. 



— + 



or 



-■-^ + \\p,\\ (1 - \\p,\\) + Cie^" + C2e^\\oge\ > 

WPef - I \\Pe\\ + ^ - Ci£" - C2£2-"|log£| < 





which implies that, for e sufficiently small, 

||PelLi(x,.)>l-Co,e"-a,,£2-"|loge|. 
The result therefore follows from the normalization of in L^{Bi). 

If a > 2, the upper bound contained in the proof of Theorem l2.2l gives immediately the following bound 

dfr^ > l-Ce^llogel 



and then, using the normalization of the result is proved. 



□ 



Proof of Proposition 12.51 

Let us first consider the case < a < 2: as in the proof of Eq. (|2.15p . we first need to prove a pointwise 
estimate for = IvP^^I and, in order to find such an upper bound, we have to estimate the chemical 
potential, appearing in the variational equation satisfied by '^f^. From the definition of the chemical 
potential and the upper bound contained in the proof of Theorem 12.21 we immediately get 

e Pe < E^ +Cie +62£|loge|+e ||w^ \\L-i{Bi) 
but, from the same upper bound we obtain 



and then 



< 



Cie2" + C2e'|loge| 



4 

By replacing the above bound in the variational equation, we get 

16e" 



(3.40) 



1 02 



Cie^" + Gael log £|-C3£'"{/, 



and we can conclude that there exists a constant c (depending on fJi), such that the function C/e(r) is 
subharmonic for any < 1 — cs" . Following again the proof of Eq. (|2.15p , we can therefore obtain the 
following estimate 

^ C\\Ue\\LHBA-D.) 

for any r E Bi such that 

r < \/l - ce° - g. 

Choosing for instance g = e" /3 and using (|2.24p . we can conclude that there exists a constant Cq^ such 
that 

C/e(r) <Coi epilog £| 

for any r € 7^'. The result then follows from the application of the comparison principle to the variational 
equation satisfied in 7^' by 

Cn^e 3 I loge| 

The case a > 2 can be treated in a similar way. The only difference is in the upper bound for the chemical 
potential (|3.40p . which in this case becomes 

The subharmonicity of can now be proved in the region r < \ — ^ for any 1 < /3 < 2. As a 
straightforward consequence of (|2.25p , we then get the pointwise estimate 

C/e(r)<Cn,e^|loge| 
for any r G 7^" . The result is again obtained by means of the comparison principle. 



□ 



4 Conclusions and Perspectives 

We have analyzed rigorously the leading order asymptotics for the ground state energy and the density 
profile of a rapidly rotating Bose-Einstein condensate in a fiat trap with a finite radius in the limit where 
the coupling parameter is large. Depending on the scaling of the rotational velocity with the coupling 
parameter, different asymptotic density functionals emerge. 

Our estimates are based on trial functions that capture the essential features of the expected vortex 
structure and show the possible formation of "holes" where the density is exponentially small as a function 
of the inverse coupling parameter. The error terms in our estimates are of the expected order but the 
bounds are not sharp enough to exhibit the details of the fine vortex structure. Nevertheless, we can 
prove that rotational symmetry is broken in the ground state for const.] loge] < r2(e) < const. /e. 



An important open problem is to carry the analysis further to the next to leading order and investigate 
the transition of the vortex lattice to a "giant vortex" at high rotational velocities. 
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